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Abstract. Deep Inelastic Scattering and Drell-Yan experiments have measured a light flavor asymmetry in the proton sea. 
The excess of d over u quarks can be understood in many models, but the ratio d(x)/u(x) measured by Fermilab E866 has 
not been successfully described. Fermilab E-906 will probe the kinematic dependence of this ratio with better resolution and 
extend it to higher x. We have developed a hybrid model that includes both perturbative and non-perturbative contributions to 
the proton sea. A meson cloud formalism is used to represent the non-perturbative fluctuation of the proton into meson-baryon 
states. We include perturbative processes by using a statistical model that uses Fock states of quarks, antiquarks and gluons to 
represent the parton distributions of the 'bare' hadrons in the meson cloud. We compare our results to the E866 data. 
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INTRODUCTION 

Asymmetry in the proton sea has been well-established by experiment [1, 2, 3]. Many theoretical models for the 
asymmetry have been proposed, based on Pauli-blocking, a meson cloud, chiral perturbation theory, instantons, or 
statistical processes. For reviews, see Kumano [4], and Peng and Garvey [5]. In this paper we examine both perturbative 
and non-perturbative contributions to the proton sea, and show that either will lead to an asymmetry. Since both 
contributions are present, we propose a hybrid model, and compare to experiment. 



NON-PERTURBATIVE CONTRIBUTION: THE MESON CLOUD 

As first suggested by Thomas [6], and later by Henley and Miller [7], the meson cloud of the proton provides a natural 
explanation of the asymmetry. The fluctuation of the proton (uud) into a neutron (add) and a n + (ud) creates an excess 
of d over u. The meson cloud is represented in a Fock state expansion of the proton in terms of mesons and baryons 

\p) = Vz\ /?}bare + £ / dy d 2 k ± fo M (y,k\) | B(y,k ± )M(l -y,-k ± )). (1) 

MB J 

The leading term in the expansion is the 'bare' proton, taken to consist of valence quarks alone or valence quarks 
plus a symmetric light sea. Z is a normalization constant, <J>bm {y, k\ ) is the probability amplitude for finding a physical 
nucleon in a state consisting of a baryon B with longitudinal momentum fraction y and meson M of momentum fraction 
(1 — y), each with squared transverse relative momentum k\. The quark distribution functions q(x) are given by 

q{x) = qbare(x) + 8q(x) , (2) 

with #bare(*) the quark distribution function of the 'bare' proton, and 

5q{x)= LU lfMB{y)qM ^)j + IJ BMiy)qB ^) d j) ■ (3) 

The splitting functions fMB(y) depend on coupling constants, form factors, and cutoffs. Models differ in the number 
of meson-baryon states included in the expansion, the parton distributions qM and qs, form factors, cutoffs and 
coupling constants. The dominant terms in the expansion are p -> nn + and p — > p7t°. The first term causes a sea 
flavor asymmetry. The second term contributes equal numbers of d and u, so it will affect the ratio dj u, but not the 
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FIGURE 1. Meson cloud calculation for the proton sea asymmetry. Two cutoffs are used, Amm = 0.92 for nucleon-meson terms, 
and = 0.98 for Delta-meson terms. For these cutoffs, the 'bare' proton probability Z = 0.58, < Nk >= 0.20, < A7t >= 0.11, 
and higher order terms (Np,Ap,pa,pa,prj) represent 11% of the proton. The data is from E866 [3]. 



difference (d — u). We have shown [8] that the inclusion of the (0 isoscalar contribution improves agreement with 
experiment, and Huang et al. [9] find that an improved description of the data at higher x is obtained if the a is 
included as well. However, the theoretical ratio still remains above 1 for large x. A full meson cloud calculation should 
include all higher order terms in the expansion, p — > An,Np,Ap,p(0 7 ... In Fig. 1 we show our calculation with only 
2 cutoff values, determined by a fit to the experimental integrated asymmetry [3] of 0.1 18 ±0.012 . We can improve 
agreement with experiment by further adjustment of cutoffs, but the ratio remains above 1 for larger values of x. 



PERTURBATIVE CONTRIBUTION: STATISTICAL MODEL 

Zhang et al. [10] have used a simple statistical model which is based on a Fock state expansion of the proton in terms 
of its quark and gluon states: 

\p >= J^c ijk \{uud}{ijk} >, (4) 

ijk 

with i the number of uu pairs, j the number of dd pairs and k the number of gluons. The probability of finding a proton 
in the state \{uud}{ijk} > is pyj = |c,^| 2 . Taking into account the processes q o qg, g o qq, and g o gg, detailed 
balance between any two states leads to 

Pijk = 2 ^ 3+2/ + 2;+« 

pooo /!(/ + 2)!j!(;>l)! f fJ )( 3 + 2/ + 2;)(n+l) + ^' 

and an integrated asymmetry of 0.123, in remarkable agreement with experiment. Parton distribution functions are 
determined by a Monte Carlo simulation of the distribution of momenta among the n partons in each Fock state. In 
Fig. 2 we show our calculation of the sea asymmetry, using this model at the starting scale and our evolution to the 
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FIGURE 2. Statistical model calculation of the proton sea asymmetry. The thin line is calculated at the starting scale of the 
model, the thick line is calculated for parton distributions evolved to 54 GeV 2 , the scale of the E866 experiment [3]. 
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FIGURE 3. Preliminary hybrid model calculation for the proton sea asymmetry. Two cutoffs are used, Ajva/ = 0.83 for nucleon- 
meson terms, and A^m = 1.05 for Delta-meson terms. For these cutoffs, the 'bare' proton probability Z = 0.63, < Nk >= 0.15, 
< An >= 0.15, and higher order terms (Ap,Np,pa,pco,prj) represent 8% of the proton. The data is from E866 [3]. 



experimental scale. As for the meson cloud model, the difference is reasonable, but the ratio is relatively constant for 
larger x, in disagreement with experiment. 



A HYBRID MODEL 

We have proposed that since both perturbative and non-perturbative processes contribute to the proton sea, a hybrid 
model should be used to incorporate both mechanisms. We use the statistical model to represent the the 'bare' proton 
parton distribution in the meson cloud model, which introduces an additional sea asymmetry. Our preliminary results 
are shown in Fig 3. Only 2 cutoffs were used, determined by a fit to the integrated asymmetry. The difference is well-fit, 
but the ratio is too flat. Calculations with further variation of the meson-baryon cutoff values are in progress [11]. 



SUMMARY 

Meson cloud, statistical and hybrid models describe the ratio d(x)/u(x) for x < 0.2, but fail to agree with the exper- 
imental trend to return to symmetry, or even u dominance at larger x. Higher statistics results from E-906/SeaQuest 
[12] should determine the magnitude of the disagreement between current theory and experiment. 



ACKNOWLEDGMENTS 



This work has been supported in part by the RUI program of the National Science Foundation, Grant No. 0855656. 



REFERENCES 

1. P. Amaudruz et al, NMC Collaboration, Phys. Rev. Lett. 66, 2712-2715 (1991). 

2. A. Baldit et al, Phys. Lett. 332, 244-250 (1994). 

3. E866 Collaboration, E. A. Hawker et al, Phys. Rev. Lett. 80, 3715-3718 (1998); J. C. Peng et al, Phys. Rev. D 58, 092004 
(1999); R. S. Towell et al, Phys. Rev. D 64, 052002 (2001). 

4. S. Kumano, Phys. Rep. 303, 183-257 (1998). 

5. J. -C. Peng and G. T. Garvey, Prog. Part. Nucl. Phys. 47, 203-243 (2001). 

6. A. W. Thomas, Phys. Lett. B 126, 97-1 1 (1983). 

7. E. M. Henley and G. A. Miller, Phys. Lett. B 251, 453-454 (1990). 

8. M. Alberg, E. M. Henley and G. A. Miller, Phys. Lett. B All, 396-399 (2000). 

9. F. Huang, R. -G. Xu and B. -Q. Ma, Phys. Lett. B 602, 67-72 (2004). 

10. Y. -J. Zhang, B. Zhang and B. -Q. Ma, Phys. Lett. B 523, 260-264 (2001); Y. -J. Zhang, B. -S. Zou and L. -M. Yang, Phys. 
Lett. B 528, 228-232 (2002); Y. -J. Zhang, W. -Z. Deng and B. -Q. Ma, Phys. Rev. D 65, 1 14005 (2002). 

11. M. Alberg and T. Matossian, in preparation. 

12. P. E. Reimer, E-906/SeaQuest Collaboration, J. Phys.: Conf. Ser. 295, 01201 1 (201 1). 



